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Considerable use has been made of the geometric-optical propagation theory to describe 
long wavelength terrestrial radio waves between the ionosphere and the ground bv Bremmer 
[1949], Wait and Murphy [1957], Wait [1961], Wait and Conda [1961], and Johler [1961, 
[1962]. The physical interpretation of pulse signals propagated around the terrestrial 
sphere as an application of the theory has been described by Johler [1962, 1963a, 19631)]. 
Indeed, the use of the Loran-C radio navigation-timing system to study the propagation of 
LF ionospheric waves has given further impetus to such a study. 

The geometric series from which the geometric-optical ray limit is derived like the full 
wave residue series or the series of zonal harmonics [Johler and Berry, 1962] is rigorous. 
However, the comparatively gross approximations usually employed in the geometric-optical 
ray limit are worthy of careful scrutiny at this time since the rigorous residue series methods 
have been worked out in some detail for large scale computers [Johler and Berry, 1962, 1963]. 

It is apparent from the rigorous theory that the geometric-optical rays do not merely 
diffract around the curvature of the terrestrial sphere but at considerable distance beyond 
the geometric-optical horizon are built up to stronger fields by additional ionospheric re- 
flections of the nature of a detached mode of propagation not ordinarily taken into account- 
by the saddle point method of the ray limit. 



1. Introduction 

The display of time-separated pulses on a low frequency Loran-C radio navigation-timing 
oscillogram [Johler, 1962] at considerable distance from the transmitter can be explained by the 
geometric-optical rays, figure 1. Thus, the first pulse to arrive at the receiver (o) from (he 
transmitter (s) can be considered as a signal propagated via the ground wave (j = 0) over tin 
surface of the ground with which the geodetic d=ad is associated (the waves in the earth are 
completely absorbed and hence do not ordinarily emerge again). Another pulse will arrive 
20 to 40 /zsec later by propagating over the ray = 1) as a single ionospheric reflection. Higher 
order (j = 2, 3, 4 ... ) pulses will arrive at the receiver (0) by propagation over rays which 
reflect two or more times between the ground and the ionosphere, arriving at the receiver at 
ever later times. These pulses are quite calculable for measured models for the ionosphere 
[Johler, 1963a and 1963b] from the Fourier transform-integral theorem, employing the theory of 
propagation for either the geometric-optical series or the more conventional residue (mode) 
series. Concerning the latter, Bremmer [1949] in his well-known textbook is quoted, "We 
begin by remarking that for the skywave the development of the residue method is, in general, 
less clear than the geometric-optical method . . . ." Indeed, the geometric-optical method 
provides a clear physical interpretation of the rays without the pulse solution, where, on the 
other hand, the residue or mode theory does not split into time-separated pulses without the 
Fourier inversion to the time domain. Thus, the terms of the residue series (modes) are not 
clear as far as a physical interpretation of the time domain is concerned. Both methods of 
course should yield identical numerical values, but the approximations used in the geometric- 
optical ray limit of the rigorous series are subject of scrutiny in this paper by a comparison 
with the rigorous residue or zonal harmonics methods. 



i This work was sponsored by the Advanced Research Projects Agency under ARPA Order No. 183-62, Amendment 4, or NBS Project 854 J 1 
initiated under this contract. 
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Figure 1. Diagrammatic illustration of the geometric 
optical rays in the ground-ionosphere guide. 

2. Theory of the Geometric-Optical Series 

The propagation problem described, figure 1, has been formulated rigorously for a source s 
and an observer o, in the guide between the ionosphere and ground without regard to geometric- 
optical rays, in a companion paper, Johler and Berry [1963]. This solution is the series of 
zonal harmonics, which for the particular case of a vertical electric field, E r , volts/meter, figure 1, 
assuming a source dipole current-moment, I Q l = Airln c^ 3. 38 (10~ 2 ) ampere-meters, 



X 
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k\ and a are constants, ki=-rii^~~ 9 a is the radius of the terrestrial sphere at a frequency/=co/27r 

c c 

where c is the speed of light. P n (z) is the solution of Legendre's differential equation, 

(1-2 2 ) d ^fi 1 -2z d ^+n(n+l)PM = 0- (2) 



Also, 



i(2)=y|j,«(2), (3) 



r, w (*)= 



-^ m% H ( Z ), (4) 



where J n+V2 (z) and H$?ly 2 (z) are Bessel and Hankel functions of order ra+K and argument (s) 
and the Hankel function H^ is of the second kind. The quantities B n and T„ are related to 
the spherical reflection coefficients, R' n and Tf,, respectively, 

t> _ —4 , n(ha) ps , . 



ff(W) T , ; (6) 



where 
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where d=a-{-h and h is the height of the ionosphere boundary above the surface and 



h 



---Jl-i^t^ (9) 



for a plasma frequency squared o? N and a collision frequency v 

The denominator in (1) is expanded in the geometric-optical series [Bremmer, 1949], 

[i-R n T n }- i =\+itRm, (io) 

provided |7?^7^|<^1, and the field E r is expressed in the geometric-optical series, 

E r =E r , +^E rJ (ll) 

i=i 

where, the ground wave, E Tt0 , which can be regarded as the zero order (j = 0) term of the series, 

1 

#r.o=TO S w(n+l)(2w + l)P n (cos 0)ff (ha^ha) {1+R s n } (12) 

and a particular ionospheric wave, j, 
^^S^+l)(2, + l)P w (cos <> {!+«}■ { ^(M 

Xrf (ha^niha) {R s n }^ { 2?} '. (13) 

The expressions (12) and (13) can be evaluated directly at long wavelengths (<30 kc/s) by 
summing the slowly converging series according to methods described by Johler and Berry 
[1962]. 

On the other hand, highly convergent representation of (12) and (13) can be obtained 
from the contour integrals, 

Er ^W L ^"^os" ^~ C ° S 6) ^h(M*.-*(*i«)[1+#-k]* (14) 

x {^^-iS!^} ,J «- 1 ^- 1 * (15) 

where the poles of order j are associated with T\- y% as the roots of v=v S) of (8), 

i ti-yAhd) , i f.»i(W). 
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Integral (14) contains the real axis poles, v=n~-=-i -> - . . ., of the original series of 

A A A A 

zonal harmonics; and the only other poles of order (j+1) are those of R s v -y 2 , or, eq (8), the roots 
v=v s of 

h f f- m (M k 2 yp v - V2 (k 2 a) ' 
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Equations (16) and (17) are readily evaluated by methods of Johler and Berry [1963] in a com- 
panion paper. The convergence of the integrals (14) and (15) is assured for n=v—~ in (10) 

by \R j v _x A T j v _y 2 \<^\, which permits interchange of order of integration and summation ]T) 

J c j 

since the series converges uniformly with the circle of convergence, \R v _y 2 T v _i A \ = l. 
Employing the approximations, 

t>.^(Aia) + (iia)*T„ (18) 

with approximations of Debye or Hankel approximation [Johler and Berry, 1963 ; Berry, 1963], 
which substitutes Hankel functions of order % for those of complex order, Bremmer [1949] 
has introduced a determination of the roots t=t s , v=v s of (7) as the roots of Riccati's dif- 
ferential equation, 

^-25 2 r+l=0 (19) 

(It 

which, incidentally, have been tabulated [Johler, Walters, and Lilley, 1959], where for vertical 
polarization, 8=5 e , 

5 e = { -ifcj***?} { (haY^ikt/kD-l }~\ (20) 

However, as has been demonstrated in previous papers, Johler and Berry [1962] (17) can be 
evaluated directly. The integral (14) then becomes the classical series of residues, upon 
closing the contour along the imaginary axis of the complex ^-plane with the right half of the 
y;-plane, 
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where R v -y 2 =r} v _ y JD v _i A . This is the well-known residue series, representing the first term or 
"pure" groundwave term of the geometric-optical series. The line integral has been found to 
be negligible by Berry [1963] in a companion paper. Computation details for this residue 
series (21) and the zonal harmonics series (12) have been detailed, Johler and Berry [1962]. 
The evaluation of the integral (21) was carried out in a companion paper by Johler and Berry 
[1963]. 

The zonal harmonics series before splitting into a geometric-optical series can also be 
transformed, 

r ~kW J e COS VT ft " H {kia) * V - * ( l ] l +Uv ~ H ] 

\ 1+ f ~" + V - H (hd) r» M (M jy *.-*(M)«Sh(*i<o v ~ h v ~*j { ] 

where the integration to be considered is performed in the complex z;-plane. 

Equation (15) represents the ionosphere waves. Indeed, the sum of these waves, 
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in which the poles of i?J_^=i; ! ,_^/Z? t _^, D p _ y „=--0 produce a residue series which is precisely the 
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negative of the groundwave. Thus, the residue series resulting from the zeros v=v s of 

again yields the total field-ground and ionospheric wave. 

It should be quite possible to evaluate (15) for each ionospheric wave, j=l, 2, 3 . . ., 
as a series of residues v=v s for the poles of B s v _y 2 =r) v _ y JD v _y 2 and T s v _y 2 =rjl_ y JDl_y 2 for D V _ H =0. 
But, these are not simple poles, hence the application of the Cauchy residue theorem becomes 
quite complicated. The first ionospheric wave is of particular interest or (23) is written as a 
series of residues for j=l, employing the integration contours used in (21), 
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where the line integral Berry [1963], 



V, 1 
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can be aeglected. Z?»_^ is identical with (17) and can be evaluated by previously described 
methods of Johler [1961] for the first series of (25). The contribution of the second series 
as a result of the poles v=v s of D l v _y 2 =Q, does not occur in Bremmer [1949] and Wait [1961] 
formulation. 

Thus, for example, Bremmer [1949] in his well-known textbook (p. 33 of his document) 
makes the substitution ^ w =if„ (using exp (iut) time function). This can be accomplished 
in an exact manner, ^»=J[?» ) +f» ) ]' This latter exact substitution in (15) results in a single 
series with only groundwave type poles of importance such as those poles in the first series 
of (25), and the effect of the ionosphere poles is thrown in essence into the excitation factors, as 
was accomplished in the less exact geometric-optical theory to be derived subsequently. 
Detailed study of these matters is reserved for future work. 

Thus, 

l #_k(M) , i tf ) -'vAhd)_ 



T) 1 — - 



*,*,_„(*,<*) ' htfl^hd) 







(26) 



for v = v s in the second series of (24). Thus, it is unnecessary to make further mathematical 
development of (25) since techniques for all functions described herein have been described 

d 2 r, d 



previously, Johler and Berry [1962 and 1963] except for the ^— 2 D v _y 2 . 



Since ■=— D v . 
ov 



has 



been worked out, a continuation of this procedure, although tedious, is tractable. The signifi- 
cance of (25) and (13) is a series of residues and a series of zonal harmonics for the ionospheric 
geometric-optical wave which, like the ground wave, is rigorous. It is therefore quite tractable 
but tedious to formulate and calculate rigorous geometric-optical ionospheric waves. This 
task is reserved for future work. An examination of the more conventional ray limit of the 
rigorous geometric-optical series will now be made. 
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The ray limit of the rigorous geometric-optical series represents an evaluation of the integral 
(15) employing approximations: 

, wa e XP {|{ 1 -0y}'-?}{ 1 -ey}-'' 

/v\ 2 I \z f /v\ 2 \ 3/2 \ I 

provided 1— I- J <<1 h -*< 1— f-j f >>1 > where the former condition applies to the low 

frequency groundwave poles and the latter to the higher frequency case (100 kc/s), z=k, r 
large. Also, 

cosvt V27rsin I 4 J 

for Im v large or, say, Im v^>2, since v is exponential order. Noting the elimination of the 
poles of the integrand due to cos vw on the real axis as a result of the approximation (30), (15) 
can be written, j=l, 2, 3 . . . or^VO, 



*J^{l+Rl- H y{Rl- H }'- l {T 9 - H }> 

k x a ) 



Ef t _ exp {JT(2j+l)} r Mv 2 ~ 
kla 4 ^2r sin 6 J -co / /v- 

x«p{-«-i«[w{ 1 -(^) , } H -* ka { 1 -(^) , }1}*- (31) 

It is common practice to evaluate this integral by the saddle point method, Wait [1961], 
Bremmer [1949], 

J"°_ A(v) exp U(v)}dv^^~ A(vj) exp {/(»,) } , (32) 

where Vj are the roots, v=Vj olf'(v) = 0, or 

or 



"^{(^V'-c^y-cav-c^y},.,; 

and 

This leads to the standard geometric-optical formula, Johler [1961], using t'j=tj—Dj/c where for 
the ground wave, j=0, to=t — r)i d/c, 

E riJ =io>dDj l Cexp (iatyG'jGfaFjCj (36) 

where 

Dj^2j[(a+h) cos <Pi,j—a cos t ; ], 
and 

F,=^F;=[Fr] 2 =[l+^_ H ] 2 =[l+B] 2 (37) 
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for transmitter (t ) and receiver (r) over identical ground, 

c7 i =[B v _„p- i [r f| _„]'=[fir i m' 



(38) 



where R and T are interpreted as the well-known Fresnel reflection coefficients for vertical 
polarization (e), 

h 



R=R e = 
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(39) 



(40) 



The effect of anisotropy and the variation of electron density and collision frequency with 
altitude (nonsharply bounded ionosphere) are accounted for by the factor Cj which involves 
the reflection coefficients R, T. In the anisotropic case [Johler and Harper, 1962; Johler, 
1962], R is replaced by R e and R m , the vertical electric and vertical magnetic reflection coeffi- 
cients, and Tis replaced by T ee , T em , T me , T mmy where T ee is a vertical electric incident wave with 
a corresponding vertical electric reflected wave, T em is the abnormal component or vertical mag- 
netic resulting from the vertical electric field at the anisotropic boundary. Similarly, T mm 
refers to the vertical magnetic incident and reflected and again T me is a corresponding abnormal 
component. Thus, for the anisotropic case, Bremmer [1949], 



where 



c 1 d'V l+A x x "I 
j j\R e dx j [_l—A 2 x—A s x 2 _\ x =o' 



■^M Hm-*- mm 

-A-2 == -t^e J- ee\ Km-*- mi 

A3 = 1\ e I\ m [ J ce I n , 



- 1 - em J mc\' 



Explicit expressions for Cj are given by Johler [1961]. Johler and Harper [1962], Jollier [1962] 
have devised computation boundaries to evaluate electron density-altitude, N(h) collision 
frequency-altitude v(h) profiles of the lower ionosphere. Also, monoeiiergetic electron-ion 
collision frequencies have been introduced into these reflection coefficients, Jollier and Harper 
[1962b]. 

G?#J=# 2 ^sin 2 Tj (41) 

for identical transmitter and receiver dipoles, C = l0~ 7 I o l/d, where d=ad, the geodetic, tj=t — 
ViDj/cirj^l) and 



<*j'= 



a-\-h 



2j sin 



2j 



sin 



(a-\-h)—a cos ^-= 
(a+h) cos — .— a 



A } . 



(42) 



(A 4 = l) 



The factor aj is corrected in the caustic region, tj=-> with the factor A^t^I, Wait [I960] 



Aj^^H^( Zl ) exp {-i[5*/12-zj]}, 



(43) 



where Zj=^k x a cos 3 r,/3 sin 2 t } . Also, in the same region, Wait and Conda [1958] have replaced 
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the factor F=(l-\-R) 2 with the solution for a plane wave wrapped around an infinitely long 
cylinder, 

F- r ^T H exp — %k x aO f \ \ v v — ' / dp, (44) 

where d'=(d—d H )/2a and d H is the length of geodetic from transmitter to horizon, r^=^ and, 



Also, 

TF 1 (p)=exp[-2^/3]V7r/3(-p)^g ) {!(-p) 3 ^}, (45) 

W»=exp [-2wi/3]JJ3t>H% {*(-p)»}, (46) 

where 

7T . .57T 

~3<arg p< y 

and 

a rg ( — p)=arg p— tt, 
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arg p^ /w =— arg p where m and 7i are integers, whereupon all functions involved are single 
valued. This can also be written as a residue series summation, 

Ft--** exp M MH ± [^I^ff^ 1 - (47) 

The use of the cylindrical geodetic in the calculation of the factor F\ r in (44) can be re- 
placed by the exact solution for a plane wave on a sphere using either the series of zonal har- 
monics or the residue series of Johler and Berry [1962, 1963] or E r =E r0 (12), (14), where E 
is the primary field 

Um^=Fy. (48) 

b->oo J2JQ 

3. Discussion 

The advantages of using geometric-optical theory to interpret pulses in the time domain 
and to introduce local changes in the reflection processes warrants an investigation to determine 
corrections or extensions of the theory to make the accuracy of the predicted fields correspond 
with the same degree as the accuracy of the residue methods. To this end the computation is 
studied. The ground wave E r>0 portion of the conventional geometric-optics is calculated by 
methods developed by Bremmer. The formulation of Bremmer's [1949] theory for the large 
scale electronic computer was presented by Johler, Kellar, and Walters [1956]. A comparison 
of this method of computation, which employs some approximations, with the rigorous solution 
originally proposed by Watson [Johler and Berry, 1963] using residue methods and also using 
the series of zonal harmonics gave identical results. 

This leads one to inquire as to how well the higher order terms j=l, 2, 3 . . .of the geo- 
metric-optical series (36), or indeed, the sum of all significant terms of the series compares with 
the rigorous solution (23) (25) (16). Employing the computation methods described in the 
companion paper by Johler and Berry [1963], figures 2, 3, 4 illustrate such a comparison at 
100 kc/s, 30 kc/s, and 10 kc/s as a function of distance for various situations. A sharply 
bounded, isotropic, model ionosphere was employed for this purpose. At 100 kc/s, the geo- 
metric-opitcal ray series clings quite accurately to the rigorous curves at short distances, <C400 
km, probably because of the dominance of the Bremmer-van der Pol ground wave term in the 
series. Indeed, as distance is increased, the accuracy is still quite remarkable in spite of the 
comparatively gross approximations in (36). The ripples caused by interference of ground 
and ionospheric waves are followed quite accurately out to approximately 1,000 km. Here 
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Figure 2. Comparison of the amplitude of the field 
(E r ) as calculated with the rigorous residue series 
with the geometric-optical ray series approximation 
(100 kc/s), including diffractive corrections. 
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Figure 3. Comparison of the amplitude of the field 
(K r ) as calculated with the rigorous residue series 
with the geometric-optical ray series approximation 
(h = 67.5 km, 30 kc/s, 10 kc/s), including diffractive 
corrections. 
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Figure 4. Comparison of the amplitude of the field 
(E r ) as calculated with the rigorous residue series 
with the geometric-optical ray series approximation 
(h = 70 km, SO kc/s, 10 kc/s), including diffractive 
corrections. 
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Figure 5. Poles of contour integral in complex 
v-plane together with excitation factor, illustrating 
quantization of the propagation. 
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also the departure is not serious since no new phenomena are observed in the rigorous theory. 
The departure does, however, become serious at 2,000 km and greater distance where the field 
decays more rapidly than is indicated by the rigorous solution. A similar situation is noted, 
figures 3 and 4 at 30 kc/s or 10 kc/s. However, at 10 kc/s it is interesting to note that the seri- 
ous departure between the geometric-optical and rigorous method is smaller even at distances, 
say, 2,500 km. This, of course, is contrary to intuition, since it is generally conceded that 
geometric-optical rays are valid at the high frequency ray limit,/— > oo . However, the compari- 
son is more subtle, since, the so-called geometric-optical ray theory used at long wavelengths, 
(36), does employ spherical corrections as described previously. This leads one to inquire 
as to the cause of the departure of (36) from the rigorous theory and suggests a possible cor- 
rection to extend the theory to greater distances. 

It is possible to question the validity of the approximations for the convergence coeffi- 
cient, «j, but the calculations have been restricted to distances for which the spherical cor- 
rections are valid. The use of the cylindrical approximation (44) of Wait and Conda [1958] 
to calculate F) tT could also be questioned. However, a comparison of this factor with the 
rigorous factor (48) indicates good agreement in the distances of concern, near (d—d H =0). 
This leads one to question the saddle point approximation. 

The complex zKplane of the contour integral (16), representing the rigorous solution to 
the problem, for which computation details have been presented by Johler and Berry [1963] 
in a companion paper is illustrated, figure 8, at a frequency of 100 kc/s. The position of the 
poles of the integrand are shown graphically, and the excitation factor for each residue is also 
given. A rather interesting change occurs in the topography of the complex ^-plane at 100 
kc/s. Thus, at lower frequencies, the poles "march" downward to the left in the fourth quad- 
rant, approaching ever closer to the Im ?;-axis with ever increasing imaginary in a quite reg- 
ular manner similar to the poles with excitation factor 232, 169, 144, 132, figure 5. However, 
at 100 kc/s there is a group of poles in the region &ia<Re v<ikid which exhibit a sharp dis- 
continuity in this otherwise regular progress, such that the Im v is comparatively small and 
the excitation factor is comparatively larger, the residues of the poles with excitation factors 
286, 232, 315, 226, etc., obviously dominate the field at short distances. But, since the Im v 
is comparatively large, these modes decrease rapidly with distances and at 1,500 or 2,000 
km no longer dominate. At this point the modes with the somewhat smaller excitation fac- 
tors, 75, 120, 3.2, 2.2, etc., dominate, since their Im v, which is of exponential order, is con- 
siderably smaller. Bremmer [1949] alludes to this quantizing of the field by referring to 
these latter modes as B-poles, where the former are A-poles. The saddle point approximation 
apparently implies the A-poles as the dominant field at the surfaces of the earth. It would 
be necessary to represent the B-poles by rays skirting the ionosphere. Whereas these rays 
never touch the ground, they fill the space and apparently make a contribution at greater 
distances, (>2,000 km). 

This suggests evaluating the second summation of (25) as an additive correction to the geo- 
metric-optics. Since this residue series involves simple poles for the first ionospheric ray, it is 
comparatively simple to evaluate by previously described methods, Johler and Berry [1963]. 
The more complicated first residue series of (25) can also be evaluated for the complete rigor- 
ous geometric-optical first ionospheric reflection. It is apparent that a similar rigorous evalu- 
ation of the higher order terms of the geometric-optical series (j=2, 3, 4 . . . ) can be made 
by similar procedures, the most serious complication being the evaluation of the residues of 
higher order poles. This is complicated but tractable. 

The poles of the second series together with excitation factors 2 as calculated by (16) are 
illustrated in table 1. The poles of the diffraction factors, R v -y 21 (17), are identical with the 
ground wave. 



2 The excitation factor of the second term-series of (25) is 

v .( tf5-i)rf%(M)*, a -H(*i<0 { 1+J5-H ) 
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Table 1. Illustrating a comparison of the topography of the v -plane for the ionosphere modes of D^_i/ 2 = with 

mode B-poles of 1 — R v -i/ 2 T v -i/2 = 
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13367. 29 


-5.66* 
-5. 74i 
-5. 78i 


75 
5 
5 


13418.17 
13391.48 
13369. 91 


-5. 66* 
-5. 71? 
-6. 60f 


22 

4 

126 



Note: — i Im vi^—i (Icia) 1 ? 3 Im to^— 15 i for ground wave poles (18) which have been tabulated by Johler, Walters, and Lilley [1959]. 

The ionosphere poles are close to the B-poles of the residue series or mode series repre- 
sentation. However, the excitation factors are different. Note the large excitation factor for 
5=1 of approximately 75. 

Since the ionosphere poles have small Im v, these are not negligible at greater distance 
since the excitation factor is sufficiently large. Indeed, these residues are of sufficient amplitude 
to compensate or correct the geometric-optical curves, figures 4, 5, where the first and second 
term of the geometric-optical series are sufficient. Obviously, the corrections for the higher 
order terms of this series are required at the greater distances O~2,000 km). It can then be 
concluded that the geometric-optical rays are not merely built up by ray diffraction around 
the terrestrial sphere but indeed are built up by certain ionospheric reflection of the nature of 
a detached mode. The detailed study of these phenomena is reserved for future work. 

4. Conclusions 

The geometric-optical ray theory, which has been corrected with the conventional spheri- 
cal corrections, is a reasonably valid theory to distances of approximately 1,500 or 2,000 km. 
At greater distances and especially at LF (100 kc/s), it is necessary to correct the field for addi- 
tional ionospheric propagation, in which case an ionosphere reflection coefficient residue cor- 
rection can be made. It is quite possible with presently available techniques to evaluate the 
first term (j=l) of the rigorous geometric-optical series, and the extension to higher order 
(j=2, 3, 4 . . .) terms of the series is tractable but complicated. In any case it is clear that 
a diffraction theory is not always adequate to describe the field from a geometric-optical point 
of view, especially at distances greater than 1,500 km. The results of this analysis suggest 
such extensions for future work. 
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